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Abstract 

We show analytically that the [0,1], [1,1] and [2,1] Fade approximants of the mean cluster 
number S{p) for site and bond percolation on general d-dimensional lattices are upper bounds 
on this quantity in any Euclidean dimension d, where p is the occupation probability. These 
results lead to certain lower bounds on the percolation threshold pc that become progressively 
tighter as d increases and asymptotically exact as d becomes large. These lower-bound estimates 
depend on the structure of the d-dimensional lattice and whether site or bond percolation is being 
considered. We obtain explicit bounds on pc for both site and bond percolation on five different 
lattices: d-dimensional generalizations of the simple-cubic, body-centered-cubic and face-centered- 
cubic Bravais lattices as well as the d-dimensional generalizations of the diamond and kagome (or 
pyrochlore) non-Bravais lattices. These analytical estimates are used to assess available simulation 
results across dimensions (up through d = 13 in some cases). It is noteworthy that the tightest lower 
bound provides reasonable estimates of pc in relatively low dimensions and becomes increasingly 
accurate as d grows. We also derive high-dimensional asymptotic expansions for pc for the ten 
percolation problems and compare them to the Bethe-lattice approximation. Finally, we remark 
on the radius of convergence of the series expansion of S in powers of p as the dimension grows. 



Electronic address: torquato@electron.princeton.edu 



^Electronic address: |yjiao@princeton.edu 



1 



I. INTRODUCTION 



There has been a longstanding interest to understand the effect of dimensionahty on the 
structure and bulk properties of models of condensed phases of matter, especially lattice 
models [l-7|. More recently, the high-dimensional behavior of interacting many-particle 
systems has received considerable attention and led to insights into low- dimensional systems. 



This includes studies of models of liquids and glasses jsl-lG 



particle configurations and their local density fluctuations 



17, 



problems 



hyperuniformity of many- 



21 



19| and their relationships to classical ground states |20[, densest sphere packings 



covering and quantizer 



22] . and Coulombic systems [23|. The preponderance of studies aimed at elucidating 



the dependence of dimensionality across all dimensions have been carried out for Ising-spin 
and lattice-percolation models; see, among the multitude of such investigations, Refs. |2|-|7. 
Virtually all of this work has been carried out on the (i- dimensional hypercubic lattice "L^ . 
The present paper is concerned with the prediction of Bernoulli nearest-neighbor site and 
bond percolation thresholds on general (i-dimensional lattices in Euclidean space M"^. 

While it is well-known that critical exponents first take on their mean-field dimension- 
independent values when (i = 6, independent of the lattice, the percolation thresholds 
generally depend on structure of the lattice and are believed to achieve their mean-field 
values only in the limit of infinite dimension [l|. Whereas thresholds are known exactly for 
only a few lattices in two dimensions 2J], there are no such exact results for > 3 for finite 
d. Thus, most studies of the determination of lattice thresholds in any finite dimension have 
relied on numerical methods or approximate theoretical techniques [g, |25|-|40|. 

It has recently been shown that the [0,1], [1,1] and [2,1] Fade approximants of the 
density- dependent mean cluster number S for prototypical (i-dimensional continuum perco- 
lation models provide lower bounds on the corresponding thresholds j41]. Specifically, these 
results apply to overlapping (Poisson distributed) hyperspheres as well as hyperparticles 
of nonspherical shapes with some specified orientational distribution function. The sharp- 
ness of these bounds showed that previous simulations for the thresholds were inaccurate 
in higher dimensions, which then led to studies that reported improved estimates for the 
thresholds of overlapping hyperspheres 42| as well as for overlapping hyperparticles with a 
variety of specific shapes [43| that apply in any dimension d. 



Using the same techniques as was employed in Ref. |4l|, we obtain analogous lower bounds 
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on pc for site and bond percolation for general dimensional lattices in R . We demon- 
strate that these general lower bounds become progressively tighter as d increases and exact 
asymptotically as d becomes large. Employing these general results, we derive explicit ex- 
pressions for lower bounds on for five distinct lattices: (i- dimensional generalizations of 
the simple-cubic, body-centered-cubic and face-centered-cubic Bravais lattices as well as the 
d-dimensional generalizations of the diamond and kagome (or pyrochlore) non-Bravais lat- 
tices. Our analytical lower-bound estimates of these ten different percolation problems are 
then employed to assess available simulation results across dimensions (up through d = 13 
in some cases). We show that the tightest lower bound provides reasonable estimates of Pc in 
relatively low dimensions and becomes increasingly accurate as d grows. Our investigation 
also sheds light on the radius of convergence of the series expansion of the mean cluster 
number S{p) in powers of the occupation probability p across dimensions. 

The rest of the paper is organized as follows: We provide fundamental definitions in Sec. 
II and derive lower bounds on the percolation threshold pc in Sec. III. In Sec. IV, we describe 
the (i-dimensional lattices that will be considered here as well as obtain series expansions of 
S{p) and asymptotic expansions of the lower bounds on pc- In Sec. V, we explicitly evaluate 
the bounds on p^ across dimensions and compare them to available simulation results. We 
close with concluding remarks and discussion in Sec. VI. 



II. DEFINITIONS AND PRELIMINARIES 



A. Bravais and Non-Bravais Lattices 



A (i-dimensional Bravais lattice in R'' is the set of points defined by integer linear com- 
binations of a set of basis vectors, i.e., each site is specified by the lattice vector 

p = riiai + n2SL2 H h Ud-iSid-i + UdSid, (1) 

where are the basis vectors for the fundamental cell, which contains just one point, and 
rii spans all the integers for i = 1,2, - ■ - d. Every Bravais lattice has a dual or reciprocal 
Bravais lattice in which the sites of the lattice are specified by the dual (reciprocal) lattice 



vector q such that q ■ p = 27rm, where m = ±1, ±2, ±3 ■ ■ ■ ; see Conway and Sloane 
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for additional details. The concept of a Bravais lattice can be naturally generalized to 



3 



include multiple points within the fundamental cell, defining a periodic crystal or non- 
Bravais lattice. Specifically, a non-Bravais lattice consists of the union of a Bravais lattice 
with one or more translates of itself; it can therefore be defined by specifying the lattice 
vectors for the Bravais lattice along with a set of translate vectors that define the basis 
(number of points per fundamental cell). 



B. Connectedness Criterion 

Consider a d-dimensional lattice A in in which each site is occupied with probability 
p in the case of site percolation or in which each bond is occupied with probability p in 
the case of bond percolation. The lattice A can either be a Bravais or non-Bravais lattice. 
We consider Bernoulli percolation with a nearest-neighbor connectivity criterion for either 
site or bond percolation on A in which the coordination number z/<^ is the number of nearest 
neighbors. The following indicator function defines this connectivity criterion: 

J 1 if sites (or bonds) i and j are occupied nearest neighbors, , . 

^^"■^^■^^1 0, otherwise 

where Fj^ is the displacement vector between sites (or bonds) i and j. In the case of site 
percolation, 

/(rij) =Zs = za (3) 
where is the coordination number for the lattice A. In the case of bond percolation, 

J2 /(ri,) = Zb = 2(2A - 1) = 2{z, - 1), (4) 
i=i 

where it is to be noted that generally Zb > Zg for any d >2. 



C. Connectedness Functions 

The mean cluster number (or mean cluster size) S is the average number of sites (bonds) 
in the cluster containing a randomly chosen occupied site (bond). The pair- connectedness 
function ^2(1") is defined such that p'^P2{r) gives the probability that a site (center of a bond) 
at the origin and a site (bond center) j located at position r are both occupied and belong 
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to the same cluster. Essam showed that the mean cluster number is related to a sum over 



the pair-connectedness function 



3| 



S = l+pJ2P2{r). (5) 

r 

This relation can be equivalently expressed in terms of the Fourier transform P(k) of P(r): 



S = l+pP(k = 0). 



(6) 



Using the Ornstein-Zernike equation 45|| that defines the direct connectedness function C(r): 



P{k) = C{k)+pC{k)P{k), 



(7) 



where C(k) is the Fourier transform of C(r), we also can express the mean cluster number 
as follows: 

S=[l-pC{0r\ (8) 

Since -P(r) becomes long-ranged (i.e., decays to zero for large r slower than 1/r'^), S diverges 
in the limit p p^ , and hence we have from ([8]) that the percolation threshold is given by 



pc = [cm-'- 

It is instructive to note that the real-space equation corresponding to relation ([7]) is 

P(ri2)=C(ri2)+p 5^C(ry)P(r2,). 

i=i 

The sum operation here is the analog of the convolution integral in Mf^. 
It is believed that S obeys the power law 



(9) 



(10) 



S (X {pc- p) ^, 



p^Pl 



(11) 



in the immediate vicinity of the percolation threshold. In this expression, 7 is a univer- 
sal exponent for a large class of lattice and continuum percolation models in dimension d, 
including not only Bernoulli lattice and spatially uncorrelated continuum models, but cor- 



related continuum systems 



29 



30 



46|. For example, 7 = 43/18 for d = 2 and 7 = 1.8 for 



(i = 3. It is believed that when > c/c = 6, where rf^ is the "critical" dimension, the lattice- 
and continuum-percolation exponents take on their dimension-independent mean-field val- 



ues, 



29 



30 



46| which means in the case of (|TT1) that 7=1. These mean-field values are 



obtainable exactly from percolation on an infinite tree, such as the Bethe lattice for which 
Fisher and Essam [l| showed that the threshold is given by 

Pc = (12) 

za-1 

The dimensionality of the Bethe lattice is effectively infinite and therefore it is generally 
assumed that pc for (periodic) lattices approach the Bethe-lattice approximation (fT2l) in the 
limit d — )■ oo. We will see in Sec. IIV CI that this assumption is generally not exactly true. 
Note that for the large class of periodic lattices in which the coordination number z\ grows 
monotonically with d, the high-dimensional Bethe approximation becomes 

{d^oo). (13) 

za 



D. Cluster Statistics 



A k-mer is a cluster that contains k sites or bonds. The cluster-size distribution is 
the average number of /c-mers per site (bond). Thus, the probability that an arbitrary site 
(bond) is part of a k-mer is kuk, and hence 

oo 

^knk=p, p < Pc- (14) 

k=l 

Since the quantity kuk/T^kkuk is the probability that the cluster to which an arbitrary 
occupied site (bond) belongs contains exactly k sites (bonds), the mean cluster number S 
can be alternatively expressed as 

oo 

y^^k'^uk 

S=^ , P<Pc- (15) 

y^^kuk 

k=l 

E. Series Expansion for Mean Cluster Number S 

As indicated in the Introduction, our ensuing analysis requires partial knowledge of the 
series expansion of the mean cluster number S{p; d) for any dimension d in powers of p: 

S{p-d) = l + Y,Sm+i{d)p''\ (16) 

m=l 
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The (i-dependent coefficients Sk(d), which account for {k + l)-mer cluster configurations 
(A; = 1, 2, 3, ■ ■ ■ ), can be obtained in a number of different ways. A common way is to first 



obtain exphcit formulas for the cluster size distribution and then emp 



p expansion of 5* and thus the coefficients Sm+i of series ffTB]) 



25^, 



36 



distribution can generally be represented by the following relation: 



oy ( IT^ to get the 
44l |. The cluster size 



(17) 



k=l 



where Qkm is the number of cluster configurations (lattice animals) with size k and perimeter 



m associated with that cluster size 



29| . The basic calculation reduces to the determination 



of Qkm- In Appendix A, we provide an algorithm that enables one to obtain the explicit 
analytical expressions for the rii, 77-2, and in arbitrary dimension for both site and 
bond percolation for various (i-dimensional lattices. 

Another procedure that has been employed to ascertain the series f ll6p is to make use of 
the Mayer-type expansion of the pair connectedness function P{r) in terms of the connectiv- 



ity function /(r) defined by (|2]) 



Ref. 



41 



45| . In order to make contact with the techniques used in 



for continuum percolation, it is useful here to map those results for the Mayer-type 



expansion of -P(r) into the appropriate results for lattice percolation. For this purpose, this 



mapping, which amounts to replacing integrals given in Ref. 



41 



with appropriate sums. 



yields the following expansion of -P(r) to first order in p for lattice percolation: 

P(ri2) = /(ri2) +P [1 - /(ri2)] 5^/(ri,)/(r2,) + 0{p'). (18) 

j 

Substitution of (fT5]l into (0) yields, after comparison to (ITB]) . the dimer coefficient as 

52(rf) = ^/(ri,) = z„, (19) 

where a = s ot b for site or bond percolation, respectively, and is related to the coordination 
number z/^ of the lattice A via either ([3]) or Similarly, the trimer coefficient are given by 



SM = $^$^[l-/(ri,)]/(ri,),/(r,, 



(20) 



where the indices j and k run through all sites (bonds). The expressions f[T9]) and f l2U]) 
for the dimer and trimer coefficients are the lattice analogs of Eqs. (24) and (25) given in 



Ref. 



41 



for continuum percolation. In Appendix B, we illustrate how to apply Eq. (120|) by 
explicitly computing 53 for site percolation on the triangular lattice in (i.e., A2). 



III. LOWER BOUNDS ON THE PERCOLATION THRESHOLD 



It has recently been shown that the [0, 1], [1, 1] and [2, 1] Pade approximant of the mean 
cluster number S, a function of the particle number density, for prototypical ci-dimensional 
continuum percolation models provide lower bounds on the corresponding thresholds 41|. 
Specifically, these results apply to overlapping (Poisson distributed) hyperspheres as well as 
hyperparticles of nonspherical shapes with some specified orientational distribution function. 
Using the same techniques as was employed in Ref. |4l|, we obtain analogous lower bounds 
on pc for site and bond percolation for general ci-dimensional lattices in M'^. 

Let us denote the [n, 1] Pade approximant of the series expansion (ITB]) of the mean cluster 
number S by Si^n,i]- This rational function for any d is given explicitly by 

Sn+2 



m=l 



5, 



m+1 



+1 J 



P 



n,l 



1 Sn+2 
1- ^ P 



for < p < 



(n) 



(21) 



where p^^^ is the pole of the [n, 1] approximant, which is given by 

Sn+l 



(n) 

Po 



for n > 0, 



(22) 



Sn+2 

and 5*0 = 1. Here we use the convention that the sum in (12T]) is zero in the single instance 

(n) 

n = 0. The claim that we make is that the pole Pq for n = 0, 1 and 2 bounds the threshold 



Pc for general ci-dimensional lattice percolation (site or bond) from below for any ci, i.e. 

Sn+l 



Pc >Po 



in) 



Sn 



for 72 = 0,1,2. 



+2 



(23) 



For the [n, 1] Pade bounds to become progressively better as n increases from to 1 and 
then to 2, it is clear that the following conditions must be obeyed: 



S^ > S, 



5*3 ^ iS'25'4. 



(24) 



A. Proof in the One-Dimensional Case 

For the one-dimensional integer lattice Z, it is trivial to show that all [n, 1] Pade approx- 
imants of S* (n = 0, 1, 2, 3, . . .) provide lower bounds on the percolation threshold. To see 
this, note the mean cluster number S in this one-dimensional case is given exactly by 

S-\^, (25) 

1 — P 

8 



and hence the percolation threshold is trivially pc = 1- Expanding this relation in powers of 
p and comparing to f|T6|) yields 

Sm = 2, for m > 2. (26) 

We see from ([22]) that 

(n)_il/2 forn = 0, 

" I 1, for n > 1. ^^^^ 

and hence these poles always bound from below or equal the actual threshold pc = 1- 

Remark. For sufficiently small d >2, all [n, 1] Fade approximants of S" (n = 0, 1, 2, 3, . . .) 
cannot be nontrivial positive upper bounds on S. For example, it is known that for d = 2, 
Sm can be negative for some sufficiently large m 



B. [0, 1] Fade Bounds 

We will begin by proving that the [0, 1] Fade approximant of the mean cluster number, 

S ^ g[o,i] = 1 ^ = ^ir^ for ^<P<z~\ (28) 
1 - 1 - — 

provides the following rigorous lower bound on the percolation threshold Pc for all d: 

Vc>vT = -, (29) 

where we have used the identity 5*2 = [cf. ( !T9|) ] and z^ is given by 2a [cf- ([3])] and 2{zf^ — l) 
[cf. (jl])] for site and bond percolation, respectively. It follows that in the high-c? limit, the 
pole pf^ for site percolation is twice that for bond percolation on some d-dimensional lattice, 
as reflected in the asymptotic expansions given in Sec. IIVCI for speciflc lattices. 



41 



usmg 



Here we follow the analogous proof given for continuum percolation given in Ref. 
the aforementioned mapping between the continuum and lattice problem. In particular, 
bounds (100) and (101) for the pair connectedness function P(r) given in that paper become 
for lattice percolation 

P(ri2) > /(ri2), (30) 
P(ri2) < /(ri2) +p[l- /(ri2)] J] /(ri,)P(r2,). (31) 



9 



Note the similarity of the lower bound (13T|) to the \ow-p expansion (fT8|) : except here P 
replaces / in the sum and inequality f l3T|l is valid for arbitrary p. Note that since 1 — /(r) < 1, 
we also have from (EI]), the weaker upper bound 

P(ri2) < /(ri2)+P J]/(ry)P(r2,). (32) 

j 

Summing inequality ( 132|) over site (bond) 2 and using the definition ([6]) for the mean cluster 
number S yields the following upper bound on the latter: 

S < (33) 

1- S2P 

Now since this lower bound has a pole at p = 82'^ = z~^, it immediately implies the new 
rigorous lower bound on the percolation threshold (l29l) for any d. It is important to note 
that this lower bound is valid for any (i-dimensional lattice A. 

Note that a stronger rigorous upper bound on P(r) can be obtained by using the lower 
bound (!30|) in the inequality (!3T]) . namely, 

P(ri2) < /(ri2) +P[l- f{ru)] /(ri.)/(r2,). (34) 

j 

Summing inequality ( l34l) over site 2 and use of (E]) and ( l20l) gives the following upper bound: 

S < i±M^. (35) 

1- S2P 

Although this lower bound on S is sharper than ( !33|) . it has the same pole and therefore 
does not provide a tighter upper bound on the percolation threshold than (|29I) . 

C. [1, 1] and [2, 1] Fade Bounds 

The [1, 1] Pade approximant of S, given by ( 12T|) with n = 1, is more explicitly given by 



for 0<p<p'i\ (36) 



1 + 


^3" 






p 




Zqi . 





S < = ^ 



1 -3 
1 p 



provides the following putative lower bound on the threshold Pc in all Euclidean dimensions: 

Pc>rf^ = |^, (37) 
10 



where p'^^ is the pole defined by (122|) and we have made use of the identity 5*2 



Aizenman and Newman 



48 



used completely different methods to prove, for the special 
case of bond percolation on the hypercubic lattice Z'^, the following upper bound on S: 

1 



S < 



and hence 



1 - 2dp 
1 



Pc 



> 



2d' 



(38) 



(39) 



It is instructive to compare these bounds (that apply only for Z"^) to the [1,1] estimates. 
Using the fact that S2{d) = Zb = 2{2d — 1) and Ssld) = 2{2d — 1)^ for bond percolation on 
the hypercubic lattice (see results of Sec. HV|) . the [1,1] estimates and ^S7\ reduce to 

1 



S < 



1 - 

Pc > 



(2rf- l)p' 
1 



(40) 



(41) 



2d - 1 

It is seen that the [1,1] estimates f l40p and (1411) for the special case of bond percolation on 
Z'' provide sharper bounds than (1381) and (15^ in any finite dimension, and tend to the same 
asymptotic bound in the limit d ^ oo. 

Similarly, the [2, 1] Pade approximant of the mean cluster number S, given by fl2Tl) with 
n = 2, is more explicitly given by 



1 + 



P 



p 



Si 

1 - -^p 

J3 



for < p < Pq 



(2) 



(42) 



provides the following putative lower bound on the percolation threshold p^ in all d: 



Pc > Po 



(2) 



^3 



(43) 



where p'^^ is the pole defined by fl22p . Since the expansion of upper bound (1421) in powers 
of p is exact through order p^, we deduce, after comparison to the exact expansion f[T^ . the 
following upper bound on the fifth-order coefficient 5*5 (c?) for any d- dimensional lattice A: 

Slid) 



S,{d) < 



(44) 



With considerably extra effort, one can rigorously prove that (137|) and (H3|) are indeed 
lower bounds on the threshold Pc- However, this is beyond the scope of the present paper. 
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and will be reserved for a future work. Nonetheless, it is noteworthy that high- dimensional 
asymptotic expansions of fl37|) and fH3l) for both site and bond percolation on the hypercubic 
lattice Z"^ provide lower bounds on the corresponding exact asymptotic expansions, as ex- 
plicitly shown in Sec. II VC II Moreover, in Sec. |Vl we will see that available high-precision 
numerical estimates of for different lattices across dimensions support the proposition 
that (137|) and (l43ll are rigorous lower bounds on pc. 



D. [n,l] Pade Approximant 

We expect that higher-order [n, 1] Pade approximants {n > 3) of S* also provide lower 
bounds on pc ioi d > 2 for n > 3 and relatively low d provided that certain conditions are 
met. One such necessary conditions is that successive coefficients Sn+i and Sn+2 remain 
positive. For example, we have directly verified that both S[3^i] and iS'[4_i] yield lower bounds 
hi d = 2 and d = 3 for a variety of site and bond problems on a variety of lattices 



onv 



36 



44| . However, as noted earlier, because we expect S'„ to become negative at some 
sufficiently large value of n for d = 2 and d = 3, S[n,i] cannot always yield lower bounds on 
Pc for relatively low dimensions such that d > 2. In the limit d — ?■ oo, we have shown that 
the Sn are all positive and hence it is possible that in sufficiently high but finite d, S[n,i] 
gives lower bounds on p^ for any n. The reader is referred to a related discussion in Sec. IVII 



IV. SERIES EXPANSIONS OF 5 FOR VARIOUS (i-DIMENSIONAL LATTICES 
A. Definitions of the d-dimensional Lattices of Interest 

In this work, we consider the (i-dimensional generalizations of the simple-cubic lattice or 
simply hypercubic lattice Z'^ as well as d-dimensional generalizations of the face-centered- 
cubic, body-centered-cubic, diamond and kagome lattices for d > 2. While the first three 
are Bravais lattices, the last two are non-Bravais lattices, as defined more precisely below. 
It is noteworthy that generalizations of these lattices are not unique in higher dimensions. 
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1. d- Dimensional Bravais Lattices 



The hypercubic 1^ is defined by 

Z'^ = {(xi, . . . , Xrf) : Xi G Z} for d > 1 (45) 

wliere Z is tlie set of integers (. . . — 3, —2, — 1, 0, 1, 2, 3 . . .) and Xi, . . . , denote the compo- 
nents of a lattice vector. The coordination number of Z^ is z^d = Id. 

A (i-dimensional generalization of the face-centered-cubic lattice is the checkerboard 
lattice defined by 

-Dd = {(xi, . . . , Xd) G Z*^ : Xi + ■ ■ ■ + Xd even} for d > 3. (46) 



Its coordination number is Zd^ = '2,d{d — 1). Note that D2 is simply the square lattice in 
M?. The checkerboard lattice Dd gives the densest sphere packing for d = 3 and the densest 
known sphere packings for d = A and 5, but not for higher dimensions 



19 



2ll- 



It also 



49|. 



provides the optimal kissing-number configurations for d = 3 — 5, but not for > 6 

In order to define the generalization of the body-centered-cubic lattice that we will con- 
sider in this paper, we must first introduce another generalization of the face-centered-cubic 
lattice, namely, the root lattice Ad, which is a subset of points in Z'^"''^, i.e., 

Ad = {(xo,Xi,...,Xd) gZ'^+^ :xo + Xi + --- + Xd = 0} for d > 1. (47) 

The coordination number of is = d{d + 1). Note that = A^, but Dd and Ad are 
not the same lattices for d > 4. It is important to stress that the fundamental cell for the 
Ad lattice is a regular rhombotope, the (i- dimensional generalization of the two-dimensional 
rhombus or three-dimensional rhombohedron. 

The d- dimensional lattices Z^, and are the corresponding dual lattices of Z*^, Dd 
and Ad- While both and ^3 are the body-centered cubic lattice, they are not the same 
lattices for d > A. Indeed, has an unusual coordination structure for d > 4 in that the 
coordination number does not increase monotonically with d. By contrast, the coordination 
number of A^ is za*^ = 2{d + 1). For this reason, we choose to consider the lattice as 
a d-dimensional generalization of the body-centered-cubic lattice. The lattice vectors ej 
of can be obtained from the associated Gram matrix G = {Gij} =< ei,ej >, where 
<, > denotes the inner product of two vectors in W^. Following Conway and Sloane, we 
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set Gii = d and Gij = —1 {i ^ j). We note that A2 = A2 is the triangular lattice in M^. 
(We say that two lattices are equivalent or similar if one becomes identical to the other 
by possibly rotation, reflection, and change of scale, for which we use the symbol =^ The 
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m. We 



lattice A"^ provides the best known covering of M in dimensions 1-5 and 10-18 
note that while A^ apparently minimizes large-scale density fluctuations (among all point 
configurations in M'^), this is not true for the corresponding problem for d = 4, where Dl is 



the best known solution 



20| 



2. d-Dimensional Non-Bravais Lattices 



The generalizations of the diamond and kagome lattices considered here were introduced 



in Ref. |18|. Specifically, since the fundamental cell for the Ad lattice is a regular rhombotope, 
the points {0} U {aj} (j = 1, . . . ,d), where denotes a lattice vector of A^, are situated 
at the vertices of a regular (i-dimensional simplex. The (i-dimensional diamond lattice Dia^ 
can be obtained by including in the fundamental cell the centroid of this simplex, i.e., 

which leads to a lattice with two basis points per fundamental cell. By construction, the 
number of nearest neighbors to each point in Dia^ is zdis,^ = d + 1, corresponding to one 
neighbor for each vertex of a regular d-simplex (d- dimensional generalization of the tetrahe- 
dron). Note that Dia2 is the usual honeycomb lattice, in which each point is at the vertex 
of a regular hexagon. 

Similar to the construction of the rf- dimensional diamond lattice, the rf-dimensional 



kagome lattice Kag^ can be obtained by p 
nearest-neighbor bond in the A^, lattice 
these lattice points are located at 



acing lattice points at the midpoints of each 
|. With respect to the underlying A^ lattice, 



"° = (49) 
Xj = u + Pj/2 

where Pj = aj — v. By translating the fundamental cell such that the origin is at xq, we can 
also represent Kag^ as Ad © {vj}, where \j = aj/2 {j = 1, . . . d). Kag^ has d+1 basis points 
per fundamental cell, growing linearly with dimension. Each lattice site is at the vertex 
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of a regular simplex obtained by connecting all nearest neighbors in the lattice, implying 
that each point possesses 2d nearest neighbors in M*^, i..e, ^^Kag^ = 2d. We note that our 
dimensional kagome lattice is equivalent to the construction discussed in Ref. [33. 



B. Analytical Formulas for the Coefficients S2{d), Ss{d) and S^ld) 

Here we provide [using the cluster-size distribution function Uk expressions given in Ap- 
pendix A and Eq. (1151) ] explicit analytical formulas for the d-dimensional coefficients 5*2 ((i), 
5*3 (d) and 5*4 (rf) associated with the series expansion of S in powers of p [cf. f lTB]) ] for general 
dimension d in the cases of the Z'^, Dd, A^, Dia^ and Kag^ lattices for both site and bond 
percolation. In seven out of these ten problems, such (i-dimensional expansions have hereto- 
fore not been given. These coefficients together with the general lower bounds given in Sec. 
mil give corresponding explicit lower bounds on pc for these ten percolation problems. 

For the hypercubic lattice Z'^, the series expansion of S in powers of p for site and 
percolation, through third-order in p, are given respectively by 

S = l + 2dp + 2d{2d - 1)/ + 2rf(4d2 -7d + A)p^ + C(p^), (50) 

S = 1 + 2{2d - l)p + 2d{2d - 1) V + 2(8rf^ - I2d^ + + 2)p^ + C(p^). (51) 

The results fl50|) and f l5T]) agree with earlier ones reported in Refs. 4 and I5I, respectively. 

For the d-dimensional checkerboard lattice (the generalization of the fee lattice), the 
series expansion of S for site and bond percolation are given respectively by 

S = l+2d{d-l)p+2d{d-l){2d^-6d+7)p'^+2d{d-l){4:d^-24:d^+57d^-53d+12)p^+0{p^), 

(52) 

S = l+2{2d^-2d-l)p+2{Ad'^-8d^+9)p'^+2{8d^-2Ad^+12d^-8d^+27d^+131d-218)p^+0{p'^) 

(53) 

For (our dimensional generalization of the bcc lattice), the series expansion of S for 
site and bond percolation are given respectively by 



S 



l + 6p+18p'^ + A8p^ + 0{p^), d = 2 

1 + 8^ + 56^2 + 248^3 + 0(^4), d = 3 (54) 
l + 2{d+ l)p + 2{d + l){2d + l)p2 + 2{d + l)(4ci2 + d+l) + C(p4), > 4 
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1 + lOp + 46p2 + I86p3 + d = 2 

S = l 1 + 14p + 98p2 + 650^3 + C'(/), d = 3 (55) 

1 + 2{2d + l)p + 2(2rf + 1) V + 2(8c/3 + I2rf2 + 3rf + ly + 0(^4^,^ ^ > 4 
For the rf-dimensional diamond lattice Dia^, the series expansion of S for site and bond 
percolation are given respectively by 

S = i + (d+l)p + d{d+ l)p^ + d\d+ l)p^ + Oip"), (56) 

S = l + 2dp + 2dY + 2d^p^ + 0{p^). (57) 

For the rf-dimensional kagome lattice Kag^, the series expansion of S for site and bond 
percolation are given respectively by 

S =l + 2dp + 2dY + 2d^p^ + 0{p^), (58) 

S = 1 + 2{2d - l)p + 2(4c/2 -hd + 2)p^ + (16rf^ - 39c/^ + 43c/ - 18)p^ + C(p^). (59) 
The expansion for site percolation agrees with the one first reported in Ref. |33 . 



TABLE I: The d-dependent coefficients Sk{d) for site percolation. For A'^, the expressions apply 
in dimensions four and higher. For all of the other lattices, the expressions apply in dimensions 
two and higher. 



Lattice 


S2{d) 


Ssid) 


Si{d) 




2d 


2d{2d - 1) 


2d{U^ -7d + A) 




2d{d - 1) 


2d{d-l){2d'^ -6d + 7) 


2d{d - l)(4(i^ - 2Ad^ + 57^2 - 53d + 12) 


A* 


2{d+l) 


2(d + l)(2d + 1) 


2{d + l){Ad^ + d+l) 


Diad 


d+l 


d{d+l) 


d^{d+l) 


Kagrf 


2d 


2d^ 


2d^ 



The d-dependent coefficients Sk{d) are also summarized in Tables HI and HTl for site and 
bond percolation respectively for various d- dimensional lattices. We note that the coefficients 
5*2 (p), S^i^p) and S4{p) for all of the (i-dimensional lattices summarized in in these tables 
satisfy the conditions (1241) and hence the [0,1], [1,1] and [2,1] lower bounds onpc progressively 
improve as the order increases. Since nearest neighbors sites in Kag^ correspond exactly 
to nearest neighbor bonds in Dia^, it is not surprising that the coefficients Sk{d) for site 
percolation on Kag^ and those for bond percolation on Dia^ are identical, as shown here. 
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TABLE II: The d-dependent coefficients Sk{d) for bond percolation. For A'^, the expressions apply 
in dimensions four and higher. For all of the other lattices, the expressions apply in dimensions 
two and higher. 



Lattice 


S2{d) 


Ssid) 


Siid) 




2(2d- 1) 


2(2^-1)2 




2{8d^ - 12^2 + 3d + 2) 




2(2^2 -2d-l) 


2(4^4 - 8d^ + 9) 


2{8d^ 


- 24d^ + 12d^ - 8d^ + 27d'^ + ISld - 218) 


A* 


2(2d + 1) 


2(2d + 1)2 




2(8^3 + 12^2 + 3d + 1) 




2d 


2d^ 




2d^ 


Kagrf 


2(2d- 1) 


2(4^2 -5(i + 2) 




(led^ - 39(^2 + 43d - 18) 



C. Exact High-d Asymptotics for the Percolation Threshold 



Here we obtain the high- dimensional asymptotic expansions of the lower bounds on 
that were obtained from the [0,1], [1,1] and [2,1] Pade approximants of S for the hypercubic 
lattice Z'^ as well as d-dimensional generalizations of the face-centered-cubic (-Dd), body- 
centered-cubic (A^), diamond (Dia^), and kagome (Kag^) lattices. While we show that eight 
of the ten asymptotic expansions agree with the high- dimensional Bethe approximation (fT3ll . 
corresponding results for site percolation on both Dia^ and Kag^ do not. 



1. d- Dimensional Bravais Lattices U^, and A*^ 

In the case of site percolation on the hypercubic lattice if, the high-dimensional asymp- 
totic expansions of the lower bounds ( l29l) . ( 1371) and ( 143|) on pc obtained from the [0,1], [1,1] 
and [2,1] Pade approximants of S are respectively given by: 

Po > ^ (60) 
This is to be compared to exact asymptotic expansion obtained by Gaunt, Sykes and Ruskin 



by to the same order: 



1 5 31 1 . 
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The tightest lower bound is exact up through order l/d? and its third-order coefficient 19/32 
bounds the exact third-order coefficient 31/32 from below, as expected. It is noteworthy 
that the leading-order term in the exact asymptotic expansion is inversely proportional to 
the coordination number z^d = Zg = 2d. This is consistent with the high-dimensional Bethe 
approximation (1131) . Moreover, the leading-order term in the asymptotic expansion, obtained 
from the [0, 1] lower bound (i.e., Pc > I/S2), always agrees with the Bethe approximation 
since S2 = z^d [c.f. Eq. f|T9l) ]. In the instance of bond percolation on the l/, the asymptotic 
expansions of the [0,1], [1,1] [2,1] Fade lower bounds respectively yield 

These results are be compared to the exact asymptotic expansion obtained by Gaunt and 
Ruskin jsj to the same order: 

1 5 19 1 / 1 \ 

Observe that the tightest lower bound in the case of bond percolation is exact up through 
order 1/d (in contrast to the corresponding site percolation bound that is exact through 
order l/(i^) and its second-order coefficient 1/4 bounds the exact second-order coefficient 
5/4 from below, as it should. As in the case of site percolation on Z^, the leading order term 
of the exact asymptotic expansion of for bond percolation on this lattice agrees with the 
Bethe approximation f[T51) [i.e., Pc ^ 1/ z^d = l/(2(i)]. 

In the case of site percolation on the checkerboard lattice D^, the asymptotic expansions 
of the obtained from the [0,1], [1,1] and [2,1] Fade lower bounds on p^ respectively yield 

1 3 11 / 1 \ 
1 3 29 / 1 \ 

These results lead to the conclusion that the asymptotic expansion of the tightest lower 



18 



bound is exact at least through order l/d? and hence 



For bond percolation on D^;, the asymptotic expansions of the lower bounds yield 

Thus, we see that these results lead to the conclusion that the asymptotic expansion of the 
tightest lower bound is exact at least through order l/c?^, implying 

Note that the exact leading order terms of the asymptotic expansions of Pc for both site and 
bond percolation on agree with the Bethe approximation f|T3l) [i.e., Pc ~ '^1 ^d^ = Vl^*^^)]- 
In the case of site percolation on A^, the asymptotic expansions of the lower bounds 
obtained from the [0,1], [1,1] and [2,1] Fade approximants of S respectively yield 

(76) 
(77) 
(78) 

These results lead to the conclusion that the asymptotic expansion of the tightest lower 
bound is exact at least through order l/rf or, more precisely. 

For bond percolation on A^, the asymptotic expansions of the lower bounds yield 

1 1 1 / 1 \ 



1 

2d ^ 


1 

2rf2 


+ o\ 




1 

2d " 


3 

Ad? 


+ o\ 


ri 


1 

2d ^ 


1 

8^2 


+ o\ 


1^3 



19 



Thus, we see that these results lead to the conclusion that the asymptotic expansion of the 
tightest lower bound is exact at least through order l/rf^, and hence 

As in all of the previous cases, we see that the exact leading order terms of the asymptotic 
expansions of Pc for both site and bond percolation on A'^ agree with the high-d Bethe 
approximation ( fT3l) [i.e., pc ~ V^a* = l/(2(i)]. 



2. d- Dimensional Non-Bravais Lattices Dia^ and Kag^ 

In the case of site percolation on the d- dimensional diamond lattice Dia^, all three lower 
bounds yield the same asymptotic expansion, 

Pc = ^+ h.o.t, (84) 

where h.o.t indicates indeterminate higher-order terms. For bond percolation on Dia^, the 
asymptotic expansions of the lower bounds obtained from the [0,1], [1,1] and [2,1] Fade 
approximants of S respectively yield 

Pc = ^,+ o(^) (85) 



2d yd"^ 

(86) 

Pc = ]+o(^\ (87) 



d \d\ 

We know the order of the correction to the leading term since this problem is identical to 
site percolation on Kag^ described below. The exact leading order terms for both site and 
percolation on Diat^ agree with the Bethe approximation (fT3|) [i.e., pc ~ iZ-^Dia^ = l/d]. 

For site percolation on the d- dimensional kagome lattice Kag^, the asymptotic expan- 
sions of the lower bounds obtained from the [0,1], [1,1] and [2,1] Fade approximants of S 
respectively yield 

Pc=^,+ o(^) (88) 



2d 

P.= i+o(Jjl (89) 
= (90) 
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We know the order of the correction to the leading term is 0{l/d'^), which we determined 



from the exact p-expansion of S through order obtained by van der Marck 



33|. In the 



case of bond percolation on Kag^, the asymptotic expansions of the three lower bounds yield 

(91) 



1 
# 

1 

1 



(92) 
(93) 



Note that these results lead to the conclusion that the asymptotic expansion of the tightest 
lower bound is exact at least through order 1/d, and hence 



Pc 



1 

2d 



O 



(94) 



While the asymptotic expansions of Pc for bond percolation on Kag^ agree with the cor- 
responding Bethe approximation [i.e., pc ~ V-^Kag^ = Vl^c?)], this is not the case for site 
percolation [i.e ., Pc ^ d ^ = l/{2d)]. The latter observation was first made by van 



der Marck 



33| . but no explanation for it was given. We will discuss this issue in Sec. VI. 



V. EVALUATION OF BOUNDS ON pc AND S, AND COMPARISON TO SIMU- 
LATION RESULTS 

Here, we explicitly evaluate the [0, 1], [1, 1] and [2, 1] lower bounds on pc [i.e., inequalities 
(l29|) . ( !37l) and ( l43l) ] for the hypercubic lattice Z"' as well as (i- dimensional generalizations 
of the face-centered-cubic (-Dd), body-centered-cubic (A^), kagome (Kag^), and diamond 
lattices (Dia^) up to dimension 13 using the results for the corresponding coefficients 5*2 (c?), 
5*3 (d) and Si{d) listed in Tables [T] and [Tll We also employ these results to ascertain the 
accuracy of previous numerical simulations, especially in high dimensions. 

In Tables UTTl and HVt we compare the lower bounds fl29|) . fl37j) and fH3l) on the percolation 
threshold pc for site and bond percolation on the hypercubic lattice Z"' up through dimension 
13 to the corresponding simulation data. It is can be clearly seen that the [n, 1] Pade bounds 
get progressively better as the order n increases. Specifically, the [2, 1] Pade provides the 
tightest lower bound on pc, which becomes asymptotically exact in the limit d — )■ oo. The 
numerical values of Pc for both site and bond percolation lie above the associated best lower 
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TABLE III: Comparison of numerical estimates of the site percolation thresholds on the hypercubic 



lattice Z to corresponding lower bounds on pc- Simulation results for d = 2, d = 3 and d = 4 — 13 
1^, 



are taken from Refs. 



38 and 



371 . respectively. 



Dimension 




from Eq. (j43]) 


p^- from Eq. (1371) 


p^ from Eq. ([29|) 


1 


1.0000000000. . . 


1.0000000000. . . 


1.0000000000. . . 


0.5000000000. . . 


2 


0.59274621 


0.5000000000. . . 


0.3333333333. . . 


0.2500000000. . . 


3 


0.3116004 


0.2631578947. . . 


0.2000000000. . . 


0.1666666666. . . 


4 


0.1968861 


0.1750000000. . . 


0.1428571429. . . 


0.1250000000. . . 


5 


0.1407966 


0.1304347826. . . 


0.1111111111. . . 


0.1000000000. . . 


6 


0.109017 


0.1037735849. . . 


0.09090909090. . . 


0.08333333333. . . 


7 


0.0889511 


0.08609271523. . . 


0.07692307692. . . 


0.07142857143. . . 


8 


0.0752101 


0.07352941176. . . 


0.06666666666. . . 


0.06250000000. . . 


9 


0.0652095 


0.06415094340. . . 


0.05882352941. . . 


0.05555555555. . . 


10 


0.0575930 


0.05688622754. . . 


0.05263157895. . . 


0.05000000000. . . 


11 


0.05158971 


0.05109489051. .. 


0.04761904762. . . 


0.04545454545. . . 


12 


0.04673099 


0.04637096774. . . 


0.04347826087. . . 


0.04166666666. . . 


13 


0.04271508 


0.04244482173. . . 


0.04000000000. . . 


0.03846153846. . . 



bound and approach the lower bound as d increases, indicating that these data are of high 
accuracy, as shown in Figure [TJ Assuming the level of accuracy claimed in the simulations, 
our tightest lower bound ( l43l) is already accurate up to three significant figures for d > 10. 

We summarize in Table |V] evaluations of the best lower bound (H3l) on the percolation 
threshold Pc for site percolation on the Bravais lattices Dd and up through li = 13 and 
compare them to corresponding simulation data when available. Observe that fHSj) already 
provides a tight bound on the numerical estimates of Pc for in relatively low dimensions 
(e.g., d = 5 and 6). Our tightest lower bound fHSl ) estimates for this lattice should provide 
sharp estimates of pc for d > 6 (where no numerical estimates are currently available), which 
become progressively better as d grows and, indeed, asymptotically exact in the high-c? limit. 
In the case of A^, only three-dimensional simulation results are available for comparison. 

The best lower bounds on p^ for bond percolation on the and lattices are compared 
to available simulation data in Table IVI[ Note that the numerical estimates of Pc for the 
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TABLE IV: Comparison of numerical estimates of the bond percolation thresholds on the hyper- 
cubic lattice Z"^ to corresponding lower bounds on pc- Simulation results for d = 3 and d = 4 — 13 
are taken from Refs. 



31 



and 



371 . respectively. 



Dimension 




from Eq. (j43]) 


p^- from Eq. (1371) 


p^ from Eq. ([29|) 


1 


1 


1 


1 


1/2 


2 


0.5000000000. . . 


0.3750000000. . . 


0.3333333333 


0.1666666666. . . 


3 


0.2488126 


0.2100840336. . . 


0.2000000000. . . 


0.10000000000. . . 


4 


0.1601314 


0.1467065868. . . 


0.1428571429. . . 


0.07142857143. . . 


5 


0.118172 


0.1129707113. . . 


0.1111111111. . . 


0.05555555556. . . 


6 


0.0942019 


0.09194528875. . . 


0.09090909091. . . 


0.04545454545. . . 


7 


0.0786752 


0.07755851308. . . 


0.07692307692. . . 


0.03846153846. . . 


8 


0.06770839 


0.06708407871. . . 


0.06666666666. . . 


0.03333333333. . . 


9 


0.05949601 


0.05911229290. . . 


0.05882352941. . . 


0.02941176471. .. 


10 


0.05309258 


0.05283957845. . . 


0.05263157895. . . 


0.02631578947. . . 


11 


0.04794969 


0.04777380565. . . 


0.04761904762. . . 


0.02380952381. .. 


12 


0.04372386 


0.04359650569. . . 


0.04347826087. . . 


0.02173913043. . . 


13 


0.04018762 


0.04009237283. . . 


0.04000000000. . . 


0.02000000000. . . 




FIG. 1: Percolation threshold pc versus dimension d for bond percolation on hypercubic lattice as 
obtained from the lower bounds (j29p . (j37p . and (j43p as well as the simulation data. 
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TABLE V: Comparison of numerical estimates of the site percolation thresholds on the checker- 



board Dd and A'^ lattices to the corresponding best lower bounds on pc- Simu 
d = 3 and d = 4 — 6 in the case of the D^, lattice are taken from Refs. |3j and 



ation results for 



Simulation results for d = 3 in the case of the A"^ or bcc lattice is taken from Ref. 



34 respectively. 
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Dimension 


FC 


D,, 

T)^ from Eq (1431) 


yc 


A*, 

from Eq (I43l) 


3 




1666666666 


2459615 


2258064516 


4 


0.0842 


07500000000 




1304347826 


5 


0.0431 


04017857143 




1037735849 


6 


0252 


02462772050 




08609271523 


7 




0.01655281135. . . 




0.07352941176. . . 


8 




0.01186579378. . . 




0.06415094340. . . 


9 




0.008914728682. . . 




0.05688622754. . . 


10 




0.006939854594. . . 




0.05109489051. . . 


11 




0.005554543799. . . 




0.04637096774. . . 


12 




0.004545825179. . . 




0.04244482173. . . 


13 




0.003788738790. . . 




0.03913043478. . . 



Dr, lattice fall slightly below the corresponding lower bound. This again illustrates the 
utility of tight bounds to assess the accuracy of numerical data of Pc- We will see that 
the lower-bound estimate of pc for both site and bond percolation on Dd converges to the 
corresponding numerical estimates most rapidly among all of the d- dimensional lattices that 
we have studied in this paper. The reasons for this behavior are presented in Sec. IVIi 

In Table IVIIl we present evaluations of the best lower bound (143!) on the percolation 
threshold Pc for site percolation on the non-Bravais lattices Dia^ and Kag^ up through 
dimension 13 and compare them to the corresponding simulation data when available. The 
results for bond percolation on these lattices are given in Table IVIIII Again, it can be seen 
that (H3l) already provides a tight bound on the numerical estimates of pc in relatively low 
dimensions (e.g., d = 5 and 6). Again, as in the cases of Dd and A*^ lattices described above, 
it is reasonable to expect that our tightest lower bound fH3|) provides sharp estimates of Pc 
for d > 7, especially in high dimensions. These results are particularly useful in the absence 
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TABLE VI: Comparison of numerical estimates of the bond percolation thresholds on the checker- 



board Dd and A'^ lattices to the corresponding best lower bounds on pc- Simu 
d = 3 and d = 4 — 5 in the case of the D^, lattice are taken from Refs. |3j and 



ation results for 



Simulation results for d = 3 in the case of the A"^ or bcc lattice is taken from Ref. 



34 respectively. 
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Dimension 


FC 


D,, 

T)^ from Eq (1431) 


yc 


A*, 

from Eq (I33l) 


3 


1201635 


09965928450 


1802875 


1467065868 


4 


0.0490 


n 04534377720 




1129707113 


5 


0.026 


02619245990 




09194528875 


6 




01715448442 




07755851308 


7 




0.01213788668. . . 




0.06708407871. .. 


8 




0.009053001692. . . 




0.05911229290. . . 


9 




0.007016561297. . . 




0.05283957845. . . 


10 




0.005600098814. . . 




0.04777380565. . . 


11 




0.004574393818. . . 




0.04359650569. . . 


12 




0.003807469357. . . 




0.04009237283. . . 


13 




0.003218849539. . . 




0.03711056811. . . 



of numerical evaluations of Pc for such higher dimensions. 

It is clear that the tightest lower bound fHSl) on is accurate enough to enable us to 
compare the relative trends of the thresholds for different lattices in any fixed dimension d. 
Figure [2] shows the best lower bound ( l43l) on pc for site and bond percolation on the five 
different d-dimensional lattices Z*^, Dd, A^, Dia^ and Kag^. For any fixed dimension, we see, 
not surprisingly, that the threshold on is minimized among all of these lattices for either 
site or bond percolation due to the fact that it possesses the largest coordination number 
Z£,^. Similarly, the local coordination structure of the other lattices explains the trends in 
their relative threshold values. Observe that in the case of site percolation, the lower bound 
on Pc for Diad is identical to that for Kag^, since the two percolation problems are exactly 
equivalent to one another (see Sec. IV. B). 

Figure [3] shows the lower bounds on the inverse of the mean cluster number as a 
function of p as obtained from the upper bounds on S (128|) . (136!) and (142!) for site percolation 
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TABLE VII: Comparison of numerical estimates of the site percolation thresholds on the Kag^^ and 
Dia^ lattices to the corresponding best lower bounds on pc, denoted by p^. Simulation results for 
d = 3 — 6 for the Kag,; lattice are taken from Ref. 



for the Dia^ lattice are taken from Refs. 



3S and 



33 . Simulation results for d = 2 and d = 3 — 6 
34 respectively. Note that in the case Kag2, 



Pc = l - 2sin(7r/18) = 0.6527036446 ... is an exact result 



2M 





Kagrf 






Dimension 


Pc 


frnm F,n 

p^ iiUlli J_jL[. y^iOJ/ 


JrC 


from Eq. ([M]) 


2 




n ^nnnnnnnnn 
u.ouuuuuuuuu. . . 


0.6970413 


0.5000000000. . . 


3 


V.OOiyO 


u.oooooooooo. . . 


0.4301 


0.3333333333. . . 


4 


n 971 '\ 




0.2978 


0.2500000000. . . 


5 


0.2084 


0.2000000000. . . 


0.2252 


0.2000000000. . . 


6 


0.1677 


0.1666666666. . . 


0.1799 


0.1666666666. .. 


7 




0.1428571429. .. 




0.1428571429. .. 


8 




0.1250000000. .. 




0.1250000000. .. 


9 




0.1111111111. .. 




0.1111111111. .. 


10 




0.1000000000. . . 




0.1000000000. . . 


11 




0.09090909090. . . 




0.09090909090. . . 


12 




0.08333333333. . . 




0.08333333333. . . 


13 




0.07692307692. . . 




0.07692307692. . . 



on the hypercubic lattice Z'^ for d = 3,8 and 13. The zero of gives the threshold and we 
also include in Fig. [3] the associated numerical estimates of Pc- These plots clearly illustrate 
that the lower bounds on become increasingly more accurate as the space dimension 
increases. This is not surprising, since all of these lower bounds become asymptotically exact 
as the space dimension becomes large. The best lower bound, as obtained from ( 112|) . gives 
a highly accurate estimate of the inverse mean cluster size already for d = 8 and essentially 
should coincide with the exact result as evidenced by the very near proximity of the zero of 
the lower bound with the numerically estimated threshold Pc- 
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FIG. 2: Percolation threshold pc versus dimension d for site and bond percolation on the d- 
dimensional lattices Z'^, D^, A^, Dia^ and Kag^ as obtained from the lower bound (|43|) . (a) Site 
percolation. Note that lower bounds for Dia^ and Kag^ are identical, (b) Bond percolation. 




FIG. 3: The lower bounds on the inverse of the mean cluster number versus p for site per- 
colation on hypercubic lattice Z'^ for (a) d = 3, (b) d = 8 and (c) d = 13 as obtained from the 
upper bound on S (j28p . ()36p and ()42p . Included in this figure are the percolation thresholds (black 



circles) obtained from the accurate numerical study of Ref. 



35. 
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TABLE VIII: Comparison of numerical estimates of the bond percolation thresholds on the Kag^ 
and Dia^j lattices to the corresponding best lower bounds on pc- Simulation results for d = 3 — 5 
for the Kag^ lattice are taken from Ref. 



34. Simulation results for d = 2 and d 



6 for 



the Dia^ lattice are taken from Refs. 



40 and 



33l . respectively. Note that in the case Dia2, Pc 



1 - 2sin(7r/18) = 0.6527036446 ... is an exact result 



24|. 



Dimension 


JrC 


Kagj 
p^ from Eq. (03]) 




Dia^i 
p^ from Eq. (gS]) 


2 


0.524404978 


0.4000000000. . . 


0.6527036446 


0.5000000000. . . 


3 


0.2709 


0.2395833333. . . 


0.3895 


0.3333333333. . . 


4 


0.177 


0.1660649819. . . 


0.2715 


0.2500000000. . . 


5 


0.130 


0.1260229133. . . 


0.2084 


0.2000000000. . . 


6 




0.1012216405. .. 


0.1677 


0.1666666666. .. 


7 




0.08445595855. . . 




0.1428571429. .. 


8 




0.07240119562. . . 




0.1250000000. .. 


9 




0.06333107956. . . 




0.1111111111. .. 


10 




0.05626598465. . . 




0.1000000000. .. 


11 




0.05061061531. . . 




0.09090909090. . . 


12 




0.04598313360. . . 




0.08333333333. . . 


13 




0.04212768882. . . 




0.07692307692. . . 



VI. CONCLUSIONS AND DISCUSSION 



We have shown that [0, 1], [1, 1] and [2, 1] Pade approximants of the mean cluster number 
S for site and bond percolation on general (i- dimensional lattices are upper bounds on this 
quantity in any Euclidean dimension d. These results immediately lead to lower bounds on 
the threshold Pc- We obtain explicit bounds on for several types of lattices: rf- dimensional 
generalizations of the simple-cubic, body-centered-cubic and face-centered-cubic Bravais lat- 
tices as well as the c?- dimensional generalizations of the diamond and kagome (or pyrochlore) 
non-Bravais lattices. We have calculated the lower bounds for these lattice lattices and com- 
pared them to the available numerical estimates of pc- The lower bounds on pc obtained 
from [1, 1] and [2, 1] Pade approximants become asymptotically exact in the high-c? limit. 
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The best lower bound, obtained from the [2, 1] Pade approximant, is relatively tight for 
3 < d < 5 and generally provides excellent estimates of pc for d > Q. While the [0,1] esti- 
mate of Pc was proven to be a lower bound here, rigorous proofs of that the [1,1] and [2,1] 
estimates are indeed lower bounds will be reserved for a future publication. However, we 
have presented very strong evidence that the latter are indeed lower bounds for the class of 
d-dimensional lattices considered in this paper. 

We have seen in Sec. V that the estimate of pc obtained from the best lower bound f H3|) 
for both site and bond percolation on Dd converges to the corresponding numerical estimates 
in relatively low dimensions most rapidly among all of the five ci-dimensional lattices that 
we have studied in this paper. This is due to the highly connected nature of D^, e.g., it 
possesses the largest coordination number z^i^ among all of the five lattices studied here. 
In addition, we have shown that the asymptotic expansions of the lower-bound estimates 
are exact through at least 1/d^ and l/c?^, respectively for site and bond percolation on D^, 
and therefore more accurate than the corresponding asymptotic expansions for the other 
lattices. This observation is consistent with the principle that high- dimensional results 
encode information about percolation behavior in low dimensions, as is also the case in 



continuum percolation 



4l|-|43l. 



Among all of the ten percolation problems that we considered in the paper, the only case 
in which the high-ci limit of the threshold Pc does not agree with the corresponding Bethe 
approximation (fT2l) is for site percolation on the d- dimensional kagome lattice Kag^. The 
usual arguments explaining the tendency of a lattice to behave like an infinite Bethe tree l| 
apply in all of the other nine cases. For example, consider bond percolation on Dia^, which 
gives Pc ~ l/d (i.e., the Bethe approximation). This is the only specific instance in which 
a bond percolation problem can be exactly mapped to a site percolation problem, namely, 
that on the kagome lattice Kag^. Therefore, while the coordination number of the latter 
-^Kagd = 2(i, the threshold pc for site percolation on Kag^ must, in any dimension, agree with 
that for bond percolation on Dia^ and hence Pc must tend to l/d [not l/^^Kag^ = ^/{'^d)] iii 
the high-(i limit. 

It was once hypothesized that the percolation threshold of a lattice corresponded to the 
radius of convergence of the series expansion for S 25|]. This hypothesis rested on the 
assumption that S had no singularities on the positive real axis for p less than the critical 
value, i.e., the coefficients 5*2, 5*3, . . . were all positive. It was shown that at sufficiently high 
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order (e.g., 19th-order), the coefficients are sometimes negative for d = 2. This imphes 
that the critical concentration does not correspond to the radius of convergence of the series 
expansion for S for d = 2, strongly suggesting that there is a closer singularity on the 
negative real axis 



In analogy with the continuum percolation results of Ref. |41|, our present results offer 
evidence that in sufficiently high dimensions, the radius of convergence of ( !T6|) for Bernoulli 
lattice percolation corresponds to pc- The fact that the putative lower bound on pc [c.f. Eq. 
( 137|) ] obtained from the [1, 1] Pade approximant of S{p) [c.f., Eq. (136|) ] is asymptotically 
exact through second order terms implies that Eq. is also asymptotically exact, i.e., 

S{p) ~ —^5—, d^oo, (95) 

with critical exponent 7 = 1 (cf. (fTTl) ]. as expected. This in turn implies that the radius 
of convergence in the high- dimensional limit corresponds to the percolation threshold pc = 
S2/S3 because all of the coefficients of of the resulting expansion of S{p) [cf. Eq. f|T6l) ] are 
all positive. Recall that for d = 1, is given by (123]) . and hence all of the coefficients Sm 
are positive. Thus, it appears that the closest singularities for the occupation probability p 
expansion of S{p) shifts from the positive real axis to the negative real axis in going from one 
to two dimensions, remain on the negative real axis for sufficiently low dimensions d > 3, 
and eventually move back to the positive real axis for sufficiently large d. 
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Appendix A: Analytical Determination of Cluster Statistics for d-Dimensional Lat- 
tices 

In this appendix, we describe the algorithm that we have used to obtain analytical ex- 
pressions for the coefficients S2{d), S3{d) and S4{d) in the series expansion of the mean 
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cluster number S in powers of the site (bond) occupation probability p for any lattice in 
high dimensions presented in Sec. IIVBI As discussed in Sec. II, S can expressed in terms 
of the cluster-size distribution function uj. [c.f. Eq. fllSp ]. Therefore, it is sufficient for us 
to determine the expressions of [c.f. Eq. f[T7|) ]. from which the series expansion of S 
can be obtained in any specific d. The general c?- dimensional coefficient Sk{d) can then be 
determined using the fact that it is a polynomial in d, i.e.. 



The coefficients k„ are determined by solving a set of linear equations in the first several 
dimensions (e.g., 2 < d < 5) such that they satisfy the explicitly known forms for 5*^ in 
these relatively low dimensions. 

Our algorithm enables us to obtain analytically the polynomials Uk by directly enumerat- 
ing all of the distinct k-mei configurations associated with a site (bond) located at, without 
loss of any generality, some chosen origin. We note that two k-mei configurations are dis- 
tinct if they contain one or more distinct sites (bonds); see Fig. IHfor simple examples. To 
the best of our knowledge, such an algorithm has not been applied before to obtain explicit 
expressions for the n^s. Our algorithm works as follows: For a given ci- dimensional lattice, 
the vectors connecting a site (bond) to all of its nearest neighbors are determined. All of the 
k-mer configurations associated with a selected site (bond) are then generated. Specifically, 
a k-mei configuration is generated from a. {k — l)-mer configuration [k > 2) by adding a 
site (bond) that is a nearest neighbor of one of the sites (bonds) in the {k — l)-mer con- 
figuration. The total number of k-mei configurations for a site {k — 1-mer configurations 
for a bond) generated in this way is {k — l)\z^j^ ^\ where z\ is the coordination number of 
the given lattice A. Although in principle this algorithm can be employed to obtain cluster 
statistics for arbitrary k, we are only interested in the cases where 1 < k < 4 here, but for 
any dimension d. 

The k-mer configurations are then compared to one another to obtain the set of distinct 
k-mer configurations. For site percolation, we find that the set of vector displacements 
between any two sites is sufficient to distinguish a pair of k-mer configurations. For bond 
percolation, a k-mer contains k bonds and 7 associated sites (e.g., 7 = /c + 1 is the k-mer 
does not contain closed loops). The latter is simply a 7-mer in the site context. A k-mer 
configuration containing k bonds can be mapped into a configuration of k points by placing 




(A-1) 



n=l 
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FIG. 4: (color online). Two pairs of distinct trimer configurations associated with site percolation 
on the square lattice 1? . Note that the numbers indicate cluster site labels rather than ordered 
labels of the sites of 1? . (a) Two distinct trimer configurations that cannot be mapped to one 
another by any translation or rotation, (b) Two distinct trimer configurations that can be mapped 
to one another by a simple translation. 

the points at the midpoints of any bond. Note that these midpoints are not sites of the 
given lattice, but rather a new "site" decoration of the lattice. The vector-displacement sets 
for both the 7-mer configurations of the sites and the configuration the mapped k points 
are required to distinguish two /c-mer configurations of bonds. In particular, a displacement 
matrix M"^ is used to distinguish a pair of /c-mer configurations, a and /3. The components 
of the matrix M'l^ are the vector displacements between two sites (points) % and j, one in 
each fc-mer configuration (point configuration). Two /c-mer configurations are identical if 
every row of M"^ has at least one component that is a zero vector. 

Figure H] shows two simple examples of how the vector-displacement matrix M"'' can be 
applied to distinguish a pair of trimer configurations (i.e., clusters of three sites) for site 
percolation on the square lattice 7? . Figure IH^a) shows two distinct trimer configurations 
that cannot be mapped to one another by any translation or rotation. The associated 
displacement matrix is given by 

(0,0) (-1,0) (1,0) 
M"^= (1,0) (0,0) (2,0) , (A-2) 
(0,-1) (-1,-1) (1,-1) 
which does not satisfy the condition that every row has at least one zero vector. Note that 
we have set the distance between two nearest neighbor sites to be unity and the entry M"/ 
is always zero since it is associated with the common origin for any /c-mer configuration. 
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Figure 111(b) shows two distinct trimer configurations that can be mapped to one another a 
simple translation. The associated displacement matrix is given by 



M 



■a/3 _ 



(0,0) (1,0) (1,1) 
(1,0) (2,0) (2,0) 
(0,-1) (1,-1) (1,0) 



(A-3) 



While the matrix does not have zero vectors in every row, the vector (1,0) is contained in 
every row, which is the translation vector that maps the two trimer configurations to one 
another. It is clear that if the translation vector is a zero vector, the two trimer configurations 
are then identical. 

Finally, for each distinct k-mei configuration, the number of vacate sites (bonds) that 
are nearest neighbors of the sites (bonds) in the k-mei configuration is determined, which 
gives the value of the associated m (i.e., the exponent associated with 1 —p term in Eq. 
f lTTj) . Since distinct k-mei configurations that can be obtained from one another by simple 
rotation or translation have the same vacancy configuration, they contribute identical terms 
to the polynomials for n^. The total number of such /c-mers gives the value of the associated 
coefficient gkm- 

For five of the ten percolation problems considered in this paper, the expressions for the 
n/c's can be explicitly written as a function of dimensionality d, which are provided here. 
Explicit expressions for Ui, n2, and 71,4 in dimensions 2 to 5 for all of the ten percolation 
problems are provided in the Supplemental Material. 

For site percolation on hypercubic lattice Z'^, the n^'s are given by 



ni = p{l — p)^"^, 
n2 = dp\l - p)^''-\ 
ns = 2d{d - l)p3(i _ p)6'i-5 + 2d{d - l)p^{l - pf'^'^, 
n4 = |d(rf - l)(rf - 2)p\l - pf'^-^ + \d{d - l)(8rf - 7)p\l - p)8^-8+ 
+Ad{d - l)p\l - p)8^-V + dp\l - pf^-^ 



(A-4) 
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For bond percolation on hypercubic lattice Z*^, the n^s are given by 

rii = p{l — pY'^^'^, 
n2 = {2d - _p)6'i-4^ 
ns = 2d{d - ly (1 - ^ 1 (ig^2 _ + liy (1 - pfd-e^ (A-5) 

n4 = - ly (1 - pf""-^ + 16(rf - 1) V(l - py^''~^+ 
+ ^(200^3 - 552rf2 + 574d - 210^(1 - pY^d-s^ 

For site percolation on (i-dimensional diamond lattice Dia^, the n^'s are given by 

rii = p{l — pY~^^, 
n, = l{d+l)p\l-pf , 

(A-6) 



714 = |c^(3c/2 - + 4^(1 - py^-^ 
For bond percolation on (i-dimensional diamond lattice Dia^, the n^s are given by 

rii = p(l — pY"^, 

^ ^ (A-7) 
ns = ld{Ad-l)p%l-pY^-^, 

m = ^d{5d - l){5d - 2)p\l - pf^-^. 
For site percolation on (i-dimensional kagome lattice Kag^, the n^'s are given by 

Til = p{l — pY'^, 
n2 = dp^{l — pY^~^, 



(A-8) 

n4 = ^d{5d - l){5d - 2)p\l - pY 



5d-3 



Note that these expressions of n^'s are identical to those for bond percolation on Dia^. 



Appendix B: Explicit Calculation of Ss Using Eq. (|20p For Site Percolation in A2 

In this appendix, we explicitly calculate S3 using Eq. (1201) for site percolation on the 
triangular lattice (i.e., A2) as an instructive illustration of how to obtain Sk directly from 
the connectivity function /. Since the function / is only nonzero for a pair of bonds that 
are nearest neighbors of one another [see Eq. (!20|) ]. it is clear that only when site k is not 
a nearest neighbor of site 1 and when site j is a mutual nearest neighbor of sites 1 and k, 
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the product in the double sum has nonzero value (i.e., unity). This also suggests that site 
k can be at most two bonds away from site 1, otherwise it cannot share a common nearest 
neighbor with site 1. 
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FIG. 5: (color online). Three-site clusters (3-mers) of the triangular lattice that contribute to the 
coefficient (5*4. (a) The lineal configuration in which sites 1 (red or dark gray) and k (blue or light 
gray) are connected by a single common nearest neighbor j (empty circles), (b) The non-lineal 
configurations in which sites 1 (red or dark gray) and k (blue or light gray) can be connected by 
two common nearest neighbors j (empty circles). 

Figure [5] shows two configurations of sites 1 and k that contribute to Eq. (A-1). In the 
first configuration (Fig. [5|l), sites 1 and k form a straight line and can be connected by the 
common nearest neighbor j in between. Due to the symmetry of the lattice, there are 6 such 
lineal configurations, each contributing 1 to ^3. In the second configuration (Fig. [Sb), each 
pair of sites 1 and k can be connected by 2 common nearest neighbors j, which form a folded 
line. Again, due to the symmetry of the lattice, there are 12 such non-lineal configurations, 
each contributing 1 to 6*3. Thus, we have 



5*3 = 6 X 1 (lineal configurations) + 12 x 1 (non-lineal configurations) = 18. 



(A-9) 



We note that both the lineal and non-lineal configurations are 3-site clusters. One might 
initially think that a simple counting of all 3-site clusters would lead to the same result. 
Although such a counting procedure would lead to the correct result for some special cases, 
such as site percolation on the square lattice, it is generally is not valid. For example, 
the equilateral-triangle 3-site clusters do not contribute to S3 here. This naive counting 
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procedure would lead to an overestimation of 5*3. 
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